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Abstract 

The effective field theory including the dissipative effect is developed based 

on the Caldeira-Leggett theory at the classical level. After the integration 

of the small field fluctuations considered as the field radiation, the integro- 

differential field equation is given and shown to include the dissipative effects. 

In that derivation, special cares should be taken for the boundary condition 

of the integration. Application to the linear sigma model is given, and the 

decay process of the chiral condensate is calculated with it, both analytically 

in the linear approximation and numerically. With these results, we discuss 

the stability of chiral condensates within the quenched approximation. 
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I. INTRODUCTION 



Recently, there has been great interests in the collective phenomena that might be ob- 
served after the heavy ion collisions, especially in the state called the disoriented chiral 
condensate (DCC) ^g. 

The quantum chromodynamics (QCD) that describes the strong-interaction physics has 
the chiral symmetry, SUl{2) ® SUr{2)^ that is explicitly broken with the order of the pion 
mass mTT- The behavior of the vacuum state can be represented by the order-parameters 
{a) = {qq) and (vr) = {q'y^q) where g is a quark field. At zero-temperature, the effective 
potential V"eflf((cr), (vr)) has a bottom circle with the radius /tt, so that this symmetry is 
spontaneously broken and the vacuum state takes (a) ^ and (vr) = 0. 

The DCC state is defined as that on the bottom circle of Vcs but with (tt) ^ 0, and 
expected to be produced after the high-energy proton or heavy-ion collision. About the 
formation of the DCC, there has been a lot of discussions, especially concerning with the 
quenching or annealing scenarios but the present situation is still controversial. 

As has been pointed out in another important problem is the decay process of the 
DCC state. Because of the explicit breaking of the chiral symmetry, the DCC state is only 
meta-stable and considered to decay into the true vacuum with radiating many number of 
coherent pions If the life-time of the DCC is very short compared with the formation 
time, it will be very difficult to observe the DCC state through the pion signal. 

The entire process of the DCC can be formulated as the dissipative system of the collec- 
tive coordinates (the order parameters (a) and (vr)) under the back-ground pion radiation, so 
that we can analyze the DCC, especially the decay of it, with the method of the dissipative 
theory in the semiclassical dynamics. 

There has been a long history to derive the equation of motion with the radiative damping 
in the classical electrodynamics and the thermo-statistical theory. The motion of the charged 
particle with the electromagnetic radiation had been much discussed since pre-quantum 
mechanical era [^]. Here we should mention that the advanced Green function for the 
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radiation was introduced by Dirac |T0| to cancel the divergent electron's radiation mass and 



was justified by Wheeler and Feynman |]TT| . 

In the thermo-statistical theory, the dissipative dynamics has been considered with the 
Brownian motion and the Langevin equation has been one of the most important tool to 
attack the problem [P^ ]: 

Mq + ^q + V'iq) = Rit), (1) 

where q and R are the particle coordinate and the fluctuating force and 7 is the dissipative 
coefficient. However, in the Langevin theory, the dissipative coefficient is the given parameter 
and can not be fixed by itself. 

In 1983, Caldeira and Leggett proposed a method to derive the dissipation from the 
microscopic theory [0]. They start with the system-plus-reservoir model (the interacting 



system of the collective coordinates and the background freedom) and, with integrating 
out the background freedom, the dissipative equation for the collective coordinates can be 
obtained. This method was formulated in the particle dynamics, and has been applied 
for the polaron motion in a crystal (acoustic polaron), diffusion of charged interstitials in 
normally conducting metals and the dynamics of Josephson junction |T^. In the latter 



case, the reduction of the tunneling probability was predicted for the macroscopic quantum 
tunneling effect between superconductor and insulator with the Caldeira-Leggett theory and 
it was really confirmed experimentally [|14|,|15| . 



The application of the Caldeira-Leggett theory to the nuclear /hadron physics is very 
interesting: especially to the phenomena with the multi-particle production such as the 
emission of meson, photon and so on. The energy of the system (the excited nucleus or the 
fireballs, for example) is considered to decrease "dissipatively" with the particle emission 
and the Caldeira-Leggett theory can be applied to them, where the reservoir corresponds to 
the emitted particle. Another interesting application is for the DCC phenomenon. The DCC 
process is described by the non-linear sigma model, and the collective coordinates and the 
reservoir are just the condensates (a) and (vr), and the radiative pion through the formation 



and decay processes. 

In those phenomena, the dynamical variables of the system are the fields corresponding 
to the emitted particles. For the application to such a process, we have to extend the 
Caldeira-Leggett theory for the field dynamics. That extension is very interesting by itself 
because it produces the dissipative field equation. 

In this paper, we formulate the Caldeira-Leggett theory in the field dynamics, and apply 
it to the DCC process. The extension of the theory for the field theory will be given in 
the first half of this paper as generally as possible. The resultant formulation is applied for 
the non-linear sigma model, and the existence of the dissipative effects is shown analytically 
and numerically. In the final part of the paper, the decay process of the DCC is discussed 
qualitatively within the present formulation. 

II. DISSIPATIVE FIELD EQUATION WITH CALDEIRA-LEGGETT THEORY 

To illustrate the application of the Caldeira-Leggett theory to the field equation, we take 
a simple system of the order-parameter $(a;) and the background field f{x) that describes 
the particle with the mass m: 



where V"($) is the effective potential for $ and G'[<l>]/ gives the interaction between $ and 



L = i(9$)^-y($) + G(<|.)/ + i[(9/)^ 




(2) 



/. Taking variations with $ and /, we obtain the field equations from (0): 




(3a) 



(3b) 



They are easily checked to satisfy the energy-momentum conservation. 



Eq. ( pE| ) is formally solved with the Green functions G{x — y; m) that satisfies 



S'^ix - y). 



(4) 
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The fluctuation / is then written byQ 




(5) 



and, substituting it into (|3aD , we obtain the integro-differential equation for $ 



6Vmx)] 6G[<^{x)] 
5$ ^ 5$ 



/ 



d'yGix-y;m)Gmy]) = 0. 



(6) 



This equation is essentially equivalent with ( |5a| ) and (pb|). 
The Green function in @ can be written by 



(7) 



where G{k\ m) is the Fourier transform of G{x — y; m) and generally takes complex value. 
In the Caldeira-Leggett theory, the imaginary component of G{k; m) can be considered as 
the dissipative term caused by the emission of particles for the background field /. It is 
suggested by the analysis of the simple Newtonian equation including the dissipative terms: 
mx + rjx = F\ After Fourier transformation, the dissipative term r]x{t) gives that of pure 
imaginary iriujx{uj) with the spectral function x{uj) for x{t). Further, we can interpret this 
dissipation as the particle emission because ImG(A;; m) includes 6^{k'^ — m^) (on-mass shell). 
The real part ReG(fc; m) is considered to give the modification of the effective potential that 
comes from the background absorption and emission, and can be absorbed into by 
re-definition. Thus we can drop it because the parameters in V^($) is adjusted from the 
experimental data. Finally, we obtain the dissipative field equation corresponding to (1): 



d'^^{x) + 



5V[^{x)] 5G[^x)] 
5$ ^ 5$ 



/ = 



(8) 



where 





^ The homogeneous part can be dropped safely because it can be absorbed with proper adjustment 
of the boundary condition. 
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Eq. @ is still ambiguous because a variety of Green functions (advanced, retarded, 
causal and so on) can be used in it and give solutions with different boundary conditions. 
We should select the proper Green functions that shows the dissipative effect when t —>■ oo 
with real /. In the application to the linear sigma model, it will be shown that the advanced 
Green function satisfies both conditions and a proper selection. 



III. APPLICATION TO LINEAR SIGMA MODEL 

We apply the formulation given in the last chapter to the linear sigma model by Gel- 
Mann and Levy WE\ and discuss the dynamical behavior of the order parameter of the SU{2) 



chiral symmetry. 

The Lagrangian of the sigma model is 

1 



^=2 



(daY + idTtY -V{a,7t), (10) 



where a and vr = (vr^, vr^, vr'^) are the sigma and the pion fields. In this paper, we discuss 
only the phenomena related with the neutral pion so that tt in ( p!OD is a single field from now 
on. The application to the charged pion goes in the same way but we have to modify the 
Lagrangian ( ]T0| ) to include the photon degrees of freedom. The effective potential V{a, it) 
is taken to the forth order of the fields: 

V{a,7r) = ^ia' + 7r'-uY-Ha, (11) 

where the parameters A, u and H are written by 

with the sigma and pion mass m^^T^ and the pion decay constant f-,^. 

The last term in (|Tl|) breaks is the chiral symmetry explicitly and produce the finite pion 
mass. It is the simplest potential that describes the explicit symmetry breaking phenomena, 
and consistent with the low-energy theorem. 

The fields a and vr can be divided into two parts: 
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(j= {a) + 5a, vr = (vr) + 5n, (13) 

where (a) and (vr) corresponds to the order parameter of the chiral symmetry for the (dis- 
oriented) condensed state (a) = and (vr) = (^|gr75g|\I/) with the quark field 
q = {u,d). (The condensed states may be described with the coherent or the squeezed one 
as quantum states. |T^) The 5a and 5ti are the fluctuations around them, and represent the 
sigma and pi meson degrees of freedoms. The decomposition ([T3[) is essentially the same 
one in Bogoliubov theory for the weak-interacting Bose liquid. In that theory, the fluctua- 
tion represents elementary excitations (phonon/roton) that cause the dissipative effects in 
superfluids fl^ . 



Minimizing the effective potential V , we obtain the vacuum state |0) with which the 
order parameters become 

{a) = (0|a|0) = U, (tt) = (0|7r|0) = 0, (14) 

and the fluctuations 5a and 57r describes the mesons with the mass m^j and m,r each other. 

Now we consider the dynamical behaviors of the disoriented-condensed vacuum state 
|\E') 7^ |0) ((cr) 7^ and (vr) 7^ 0) as the $ in the last chapter and the fluctuations 5a and 
(Jvr as the background field /. Substituting Eq. (|T3D into (0)) obtain 

L = l [{d{a)y + {d{rr)r] - V{{a), (vr)) + H{a) 
+ d{a)d5a + din) 0571 



1 r 



+ 2 



1 



{d5a)^ - ml5a^\ + - [{0511^ - ml5r 



(15) 



where we assume that the background fields are small and describe the real sigma and 
pi meson degrees of freedom, so that we dropped 0{5a^ ,5'K^)-teim.s and the 5a'^ and ^vr^ 
terms are adjusted to be the mass term with the real pion and sigma mass mg. and m^^ as 
in (|^). In this paper, we discuss the case where the condensates are on the chiral circle: 
(cr)^ + (tt)^ = z/^, SO that (cr) and (tt) can be parameterized with the chiral angle 0: 



(cr) = I/COS0, (vr) = z/sm0. 



(16) 



Substituting Eq. (|T6D into (|^), it is obtained that 



{d4>f + Hugos i 



+ vd cos (hdda + vd sin 



+ 



(17) 



Taking the variation with 0, 5a and (Jtt, we obtain a set of Euler-Lagrange equations: 



u d 6 + Hv sin 6 — 1/ sin 5(7 + z/ cos Sn = 0, 



(18a) 



((9^ + ml)6a + cos = 0, 



(18b) 



(5^ + m^)5a + u6'^ sin = 0, 



which correspond to ( pa]) and (|3ED . 

Now we can apply the formulation developed in the last chapter to (|lH 
we obtain the dissipative field equation to ([TO): 



(9^0(x) H sin — z/ sin 1 



32 r~ 



5a + u cos (bd'^Sn = 



where and (Jvr are the similar functions as (|^). 

For G{x — y;m), we take the advanced Green function 



Gadvix -~y;m) = J 



-ik{x-y) 



:i8c) 



SH) and 



(19) 



(20) 



(27r)^ fc^ — — iesgn(k^) 

The reality of Sa and (57r can be checked easily by the direct calculation. In the next chapter, 
we will show that the advanced Green function really gives the dissipative solutions. From 
(^), we can read off the imaginary part of G{k;m): 



ilmG{k; m) = iiTsgn{k^)6'^{k'^ — m 



(21) 



so that 5(7 and Stt become 



d^6a{x) = -J^ dy^ j d^y6{x-y;m^)dl cos <f){y), (22a) 
d'6n{x) = -J^ dy'' j d'y6{x-y;m^)dlsm(f){y), (22b) 



2 Jo 

where 6{x — y; m) is the invariant deha function: 

6{x -y;m) = — ^ / d'^k6\k^ - m^)sgn{k\''''^'^-y\ (23) 
i(27r)'' J 

IV. SOLUTION OF DISSIPATIVE FIELD EQUATION 
A. Asymptotic Behavior of Solutions 

We consider the asymptotic behavior of the solutions that satisfy Eq. (|T^) and show the 
dissipative nature of it. 

The order parameter 0(x) is assumed to decrease when t — > oo, so that it can be regarded 
smaU. Then we can expand Eq. ([19D about and approximate to the first order (the hnear 
approximation). Instead of 0(t,x), we use the Fourier component: 

Mt) = / d'^t^ x)e-"^^ (24) 

for which Eq. ( [I9| ) is approximated to be 

0k(t) + k20k(t) + -— / dssinu^it - s) + kVk(s)| = 0, (25) 



with = ym^ + and a = fn/i^ ~ 1.05.0 It should be noted that no correlations exist 
among different-momentum modes in (^) in the linear approximation. 



^ This slight shift from one is resulted from the deformation of the chiral circle because of the 
explicit symmetry breaking. 
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Substituting the ansatz (f)k(t) = e "-^^"^^ into (^5]), we obtain the characteristic equation 
for the index a(k): 

a(k)^ + {2k + 1/2 + a)mla{k)^ + {k^ + (1/2 + a)K^ + a] = 0, (26) 
with K = |k|/m^. One of the solutions can be written as a(k) = p(k) + iq^n): 



Pin) = ^v'2p(«)-(2«:2 + l/2 + a), g(«) = '^^2piK) + {2k^ + 1/2 + a), (27) 



where p{k) = y/t^ + (1/2 + a;)^^ + a. It can be shown easily that p(k) takes real and 
positive value when k > 0, so that the asymptotic behavior of 0k is found to be 

Mt) ~ e'P^^')' sm[q{n)t + 6], {t oo) (28) 

with the constant phase. Summarizing the whole results, we obtain the asymptotic differ- 
ential equation to (|TU|): 

0k - k^^k + 70k + sin ~ 0. (29) 

where the dissipative coefficients are 7 = 2p(/t). The momentum dependence of p{k,) is 
shown in Fig. |I|. 

B. Numerical Results 



The fields Sa and Sw defined by { \22a\ ) and ( |22b| ) satisfy the differential equations: 



(a^ + ml)5a + 92 cos = 0, {d^ + ml)5n + d'^ sin = 0. (30) 

Hence the solutions of the original integro-differential equations are given as those of three 
differential equations, ( pJ]) and (|30D. 

In this paper, we show the numerical solutions for two cases: the uniform and the 
expanding solutions. For physical quantities, we took: 

= 92.5 MeV, M = 940 MeV, = 135 MeV, (31) 

10 



and mo- = 600 MeV was used for the mass of sigma meson. With these values, the parameters 
in (pED are fixed as A = 20.0 and v = 87.4 MeV, each other. 



1) Uniform solution. 

The solution that is uniform for the space-dependence is characterized by = (and 
correspondingly, 6a = 6a (t) and 6n = 6TT{t)). In this case, eqs. (|19|) and (|30|) are reduced 



to the ordinary differential equations, which can be easily solved. The numerical results 
are shown in Fig. ^ where the scaled time ^ = m^t has been used. Herein, the damping 
oscillation behavior is easily confirmed that is proved analytically in the last chapter. 

The dissipating behaviors can be read off in the phase diagram (Fig. |^), too, where each 
lines are the phase trajectories to these solutions; the spiral pattern around the origin show 
that they behave like the damping oscillator asymptotically. 

For quantitative check, we consider the damping rigid oscillator: 

+ 70 + sin0 = O, (32) 

with the damping coefficient 7 consistent with the dissipative coefficient p{k) in (P?]): 

7 = 2p(0) = 0.7m^. (33) 

The phase diagrams for (^) are given in Fig. ^ The trajectories in Fig. |^ are found to behave 
similar with those in Fig. |, especially in the asymptotic region (close to the origin). In the 
non-asymptotic region (far from the origin), the trajectories in Fig. ^ have the modulations 
around those of the damping rigid oscillator. They come from the non-linear dissipating 
behavior in (|T^ and (^), which is more effective in the non-asymptotic region. 

Through the comparison with the damping rigid oscillator, we can also realize the com- 
plicated behaviors at (0 = ivr, = 0) in Fig. ^. They are just the turning points of the 
rigid oscillator, and the trajectories around them are changed in chaotic manner under the 
small perturbation. 



2) Expanding solution. 
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With putting cj) = ()){t 



— X 



^) in ( [T9| ) and (|30D, we get the expanding solution in 



x-direction (uniform in other directions). Originally, this type of solution is given by Blaizot 



The numerical solution is given as the rigid line in Fig. ^ where the scaled local time 
C, = itLt^t was used. In this figure, we can find that the expanding solution is dumped faster 
than the uniform solution (shown as the dotted line in Fig. ^). To realize the expansion 
effect, we study the differential equation that Blaizot and Kryzwicki solved. In our notation, 
it becomes 



where the differentiations are for ^ = m^r. The second term in (0) is proportional to the 
time-derivative of 0, and has the effect of dissipation with the time-dependent dissipative 
coefficient 1/.^. (This effect is not a real dissipation, but a smearing of (j) brought by the 
volume expansion.) In the present case, this smearing effect has an additional effect with 
the real dissipative effect (represented by the second term in (|3^ ) asymptotically), and it 
causes the faster damping in the expanding solution. The smearing effect is found to be 
more effective in the non-asymptotic region, because the effective dissipating coefficient is 
inversely proportional with the local time r. 



We have formulated the dissipative field theory with applying the Caldeira-Leggett 
method. Explicit calculations have been done for the linear sigma model and the resultant 
field equations have been shown to have the dissipative properties with both the analytic 
and numerical ways. 

As a phenomenological application, we discuss about the disoriented chiral condensate 
that is expected to appear after the high-energy hadron collision. In the standard picture, 
the chiral symmetry is considered to be broken spontaneously at zero temperature and its 



and Kryzwicki ||T9[ with no dissipative effect. 



0" + 70' + sin = 0, 



(34) 



V. SUMMARY AND DISCUSSIONS 
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order parameters take the expectation values: (cr) 7^ and (vr) = 0. The DCC is also in 
the broken phase, but is defined to be the state where the order parameters take different 
values: (a) 7^ and (vr) 7^ 0. Because of the explicit chiral-symmetry breaking, the DCC 
state has higher energy, and it should be observed as a metastable state. 

As written in the introduction, there exist many controversies about the formation pro- 
cess of the DCC state after rebreaking of the chiral symmetry, but we concentrate on the 
decay process of it in the remaining of this paper. 

We consider the neutral DCC state with (n^) 7^ and (tt^) = 0.| The main decay 
process of this state should be the vr^-radiation, so that we can apply the above-developed 
formula with regarding vr in (^) as 7r° . The life-time tl can be estimated with the dissipative 
constant p(0) in (pT]): 

TL = l/p(0) ~ 3m-\ (35) 

If we take the quenching scenario for the DCC formation, the formation time has been 
estimated to be 

tr ~ V2m-^ ~ 0.3m-\ (36) 

It tells us that the life-time tl in ( ^5|) is ten times longer than the formation time, so that 
the neutral DCC state will be enough metastable. 

In this paper, we considered the case that the order parameters move only on the chiral 
circle (a)^ + (vr)^ = so that we could not consider the DCC formation process. The 
extension beyond the chiral circle will be given elsewhere. 



^ For the charged DCC state, photon degree of freedom is absolutely important and we have to 
extend our equations to include the dissipative effect by photon radiation. 
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FIGURES 

FIG. 1. The momentum dependence of the dissipative coefficient p{k) given by (p7|). k is 
the pion-mass scaled momentum k = |k|/m7r. The p{k,) is normahzed at k = with the value 
p{0) = 0.35m^. 



FIG. 2. A series of solutions of the integro-differential equation (19) for the space-independent 



uniform case (j){x) = (pit). The ^ is the pion-mass scaled time: ^ = mj^t and (j){t) is a chiral angle 



of the order parameter. The initial conditions are given at ^ = with # = 0. 



FIG. 3. The trajectories of the integro-differential equation (|19|) for the space-independent 
uniform case (j){x) = (t){t). The (j) and ^ are the chiral angle and the corresponding velocity with 
the scaled time ^ = m^t. The initial conditions of each trajectories are chosen as they behave 
asymptotically (||) with ^ = ^ = 0, ±1, • • • , ±5). 

FIG. 4. The trajectories of the dumped rigid oscillator (^). The dissipative coefficient 
is chosen to be consistent with the asymptotic value in the integro-differential equation ([l9|): 
7 = 2p(0) = 0.7m,r. The (j) and ^ are the chiral angle and the corresponding velocity with 
the scaled time ^ = mj^t. The initial conditions of each trajectories are chosen as in |3|. 

FIG. 5. One-dimensional expanding (scaling) solutions for the chiral angle 0(^) where ^ is a 



pion-mass scaled local time: ^ = m^rT = mT^^/t"^ — x^. The rigid line is for the dissipative case: the 
solution of the integro-differential equation (|^) , and the dashed line is for the nondissipative case: 
the original Blaizot-Kryzwicki solution of (^4|). 
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